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Abstract. We consider the first family of if (curl, fi)-conforming Nedelec finite elements on 
tetrahedral meshes. Spectral approximation (p-version) is achieved by keeping the mesh fixed and 
raising the polynomial degree p uniformly in all mesh cells. We prove that the associated subspaces of 
discretely weakly divergence free piecewise polynomial vector fields enjoy a long conjectured discrete 
compactness property as p — > oo. This permits us to conclude asymptotic spectral correctness of 
spectral Galerkin finite element approximations of Maxwell eigenvalue problems. 
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1. Introduction. Identifying spectrally correct conforming Galerkin approxi- 
mations of the Maxwell eigenvalue problem [15]: seelQ u G H (curl, fl) and u> > 0, 
such that 

curl u, curl v) i3 . n . = u? (eu,v) i2(f2) Vv G H (curl, fi) , (1.1) 

(e, fj, G (L°°(fi)) 3 ' 3 uniformly positive definite material tensors) has turned out to 
be a highly inspiring challenge in numerical analysis. Obviously, eigenfunctions of 
(|1.1[) belong to Ji(curl, fi) n if o(div £ 0, fi) and the compact embedding L 2 (fi) «-> 
ff (curl, fi) n ff (div € 0, fi) [30] relates (|1.1[) to an eigenvalue problem for a com- 
pact selfadjoint operator. However, asymptotically dense families of finite elements 
in ff (curl, fi) n if (div £ 0, fi) are not known in general. 

Let us assume that a merely ff(curl, fi)-conforming family (Wp (-^)) peN of finite 
dimensional trial and test spaces Wp(Ai) C ff (curl, fi) for (|1.1|) is employed for 
the Galerkin discretization of (|1.1|) . The corresponding discrete eigenfunctions u p G 
Wp(M), if they exist, will satisfy 

u p G X}(M) := {w p G W*(M) : K,v p ) i2(!]) = Vv p G Ker(curl) fl W$(M)} . 

(1.2) 

We cannot expect Xp(M) C if(curl, il) n Jio(div e 0, f2) and, thus, a standard 
Galerkin approximation of (jl.ip boils down to an outer approximation of the eigen- 
value problem. Good approximation properties of the finite element space no longer 
automatically translate into convergence of eigenvalues and eigenfunctions. As inves- 
tigated Caorsi, Fernandes and Rafetto in [12], an array of other requirements has 
to be met by the finite element spaces, the most prominent of which is the discrete 
compactness property [3] . 

Definition 1.1. The discrete compactness property holds for an asymptoti- 
cally dense family (Wp(A^)) p£N of finite dimensional subspaces of H (curl, fi) , if any 



*SAM, ETH Zurich, CH-8092 Zurich, hiptmair@sam.math.ethz.ch 

x We use the customary notations for Sobolev spaces like H S (Q), ff(curl, fi), H(div,Q), etc., 
and write if(curl 0, fi), £f(curl 0, fi), etc., for the kernels of differential operators. The reader is 
referred to [24, § 1.2] and [29, Sect. 2.4] for more information. 
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bounded sequence in X^{M) C i?(curl, Q) contains a subsequence that converges in 

L 2 (n). 

The same notion applies in the case of homogeneous Dirichlet boundary condi- 
tions, when (jl.ip is considered in i?o(curl, f2). In this case the eigenfunctions will 
belong to ii~o(curl, £1) n H(div e 0, f2) and zero tangential trace on 9f2 has to be im- 
posed on trial and test functions. 

The discrete compactness property of yVp(A4) is key to establishing spectral cor- 
rectness and asymptotic optimality of Galerkin approximations of (II. ip . see [11,12,15] 
for details. Small wonder, that substantial effort has been spent on proving this prop- 
erty for various asymptotically dense families of ff(curl, tt)/H o (curl, f2)-conforming 
finite elements. For the /i-version of Nedelec's edge elements Kikuchi [31-33] accom- 
plished the first proof, which was later generalized in [7,23,36], see [35, Sect. 7.3.2], [29, 
Sect. 4], and [15] for a survey. Conversely, spectral edge element schemes in 3D have 
long defied all attempts to prove their discrete compactness property, though they per- 
form well for Maxwell eigenvalue problems [15,17,39]. Partial success was reported 
for edge elements in 2D: In [9] the analysis of the discrete compactness property for 
triangular hp finite elements has been tackled, but the proof of the main result relied 
on a conjectured L 2 estimate, which had only been demonstrated numerically. The 
first fully rigorous analysis of 2D hp edge elements on rectangles was devised in [8] . 

In [29, Remark 15] an interpolation estimate was identified as crucial missing step 
in the analysis. Since then, two major advances have paved the way for closing the 
gap: 

1. In [16] M. Costabel and M. Mcintosh discovered a construction of if 1 (In- 
stable vector potentials by means of a smoothed Poincarc mapping. This will 
be reviewed in Sect. [5] of the present paper. 

2. In the breakthrough paper [21] L. Demkowicz and A. Buffa achieved a com- 
prehensive analysis of commuting projection based interpolation operators. 
To maintain the article self-contained, their approach will be explained in 
Sect. 0] and their interpolation error estimates will be presented in Sect.0 

In addition, we exploit the possibility to construct high order versions of Nedelec's 
first family of edge elements [38] by using Cartan's Poincare map [4,27-29], see Sect. [3] 
for details. Another important tool are stable polynomial preserving extension oper- 
ators developed, for example, in [1,5,22,37]. In addition, we heavily rely on spectral 
polynomial approximation estimates, see [6,37,40]. 

Thus, standing on the shoulders of giants and combining all these profound the- 
ories of numerical analysis, this article manages to give the first proof for the discrete 
compactness property of the p- version for the first family of Nedelec's edge elements 
on tetrahedral meshes of Lipschitz polyhedra f2, consult Sect. [6] for the proof. 

Theorem 1.2. The sequence (Wp(Ai)) pgN of trial spaces generated by the p- 
version of the first family of Nedelec's i?(curl, tt)- or Jfo(curl, Q)- conforming finite 
elements on a fixed tetrahedral mesh M. of a bounded Lipschitz polyhedron fl C M 3 
satisfies the discrete compactness property. 

The idea of the proof is to inspect the L 2 (S!)-orthogonal Helmholtz decomposition 
[24, § 1.3] 

w p = w p (B L 2 w", wjje Ker(curl) , (1.3) 

whose so-called solenoidal components w p belong to if (curl, ft) n i?o(div0, il). The 
above mentioned compact embedding guarantees the existence of a subsequence of 
(w p ) „ that converges in L 2 (f2). Hence, it "merely" takes to show ||w p — w^H^m) — * 
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for p — > oo in order to establish discrete compactness. Clever use of projection 
operators that enjoy a commuting diagram property, converts this task to a uniform 
interpolation estimate. The core of this paper is devoted to this seemingly humble 
program. 

Remark 1.1. Generalizations of Thm. \1.2\ to other families of tetrahedral edge 
elements, and corresponding hp-finite element schemes are straightforward [8]. For 
the sake of readability, these extensions will not be pursued in the present paper. 

Since the Poncare map does not fit a tensor product structure, extending the results 
of this paper to 3D hexahedral edge elements will take some new ideas. 

2. Poincare lifting. Let D C M 3 stand for a bounded domain that is star- 
shaped with respect to a subdomain B C D, that is, 

VaeB,xeD: {ta + (1 - t)x, < t < 1} C D . (2.1) 
Definition 2.1. The Poincare lifting R a : C°(Ti) h-> C°(JT), a e B, is defined 

as 

R a (u)(x) := [ tu(x + t(x - a)) dt x (x - a) , xeD, (2.2) 
Jo 

where x designates the cross product of two vectors in R 3 . 

This is a special case of the generalized path integral formula for differential 
forms, which is instrumental in proving the exactness of closed forms on star-shaped 
domains, the so-called "Poincare lemma", see [13, Sect. 2.13]. 

The linear mapping R a provides a right inverse of the curl-operator on divergence- 
free vectorfields, see [25, Prop. 2.1] for the simple proof, and [13, Sect. 2.13] for a 
general proof based on differential forms. 

Lemma 2.2. 7/divu = 0, then, for any a G B, curl R a u = u for all u G C 1 (il). 

Unfortunately, the mapping R a cannot be extended to a continuous mapping 
L 2 {D) i ► H^D), cf. [25, Thm. 2.1]. As discovered in the breakthrough paper [16] 
based on earlier work of Bogovskii [10], it takes a smoothed version to accomplish 
this: we introduce the smoothed Poincare lifting |f| 

R(u) := / $(a)R a (u)da, (2.3) 
Jb 

$ G C*°°(M 3 ) , supp $cB, / $(a) da = 1 . (2.4) 

Jb 



where 



The substitution 



y:=a + t(x-a) , t := y— ; (2-5) 



transforms the integral (|2.4[) into 

oo 

R(u)(x) = J J r(l - r)u(y) x (x - y)$(y + r(y - x)) drdy 

R 3 1 

= J k(x, y - x) X u(y) dy , 



(2.6) 



2 Bold symbols will generally be used to tag vector valued functions and spaces of such. 
3 The dependence of R on $ is dropped from the notation. 
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that is, R is a convolution-type integral operator with kernel 

/oc 
t(1 + r)<i>(a; + tz)z dr 

z r°° z z r°° z ( 2 ' 7 ^ 

<;<s>(x + (—)d(; + — C H(x + c—)d(. 



The kernel can be bounded by |k(aj, z)\ < K(x)\z\~ 2 , where K £ C°°(R 3 ) depends 
only on <I> and is locally uniformly bounded. As a consequence, (|2.6p exists as an 
improper integral. 

The intricate but elementary analysis of [16, Sect. 3.3] further shows, that k be- 
longs to the Hormander symbol class S'-f^R 3 ), see [41, Ch. 7]. Invoking the theory of 
pseudo-differential operators [41, Prop. 5.5] we obtain the following following conti- 
nuity result, which is a special case of [16, Cor. 3.4] 

Theorem 2.3. The mapping R can be extended to a continuous linear operator 
L 2 (D) f — * H 1 {D), which is still denoted by R. It satisfies 

curlRu = u Vu £ H(divO,D) . (2.8) 



The smoothed Poincare lifting shares this continuity property with many other 
mappings, see [29, Sect. 2.4]. Yet, it enjoys another essential feature, which is imme- 
diate from its definition (|2.2|) : R maps polynomials of degree p to other polynomials 
of degree < p + 1 . The next section will highlight the significance of this observation. 

3. Tetrahedral edge elements. In [38] Nedelec introduced a family of 
if (curl, f2)-conforming, that is, tangentially continuous, finite element spaces. On 
a tetrahedral triangulation M. of f2, the corresponding finite element spaces of degree 
p are given by 

Wp(M) := {v £ if (curl, ft): V| T £ W$(T) VT £ M} , 
Vl£(T) :={veC°°(T): v(x) = p(x) + q(x) x x, p, q G V P (R 3 ), x e T} . 

We wrote V P (M. 3 ) for the space of 3-variate polynomials of total degree < p, p £ No, 
and the bold symbol 7- > p (ffi 3 ) for vectorfields with three components in ^(K 3 ). To 
emphasize that polynomials on a tetrahedron T are being considered, we may use the 
notations V P {T)/V P (T) instead of V P {R 3 )/V{M 3 ). We also adopt the convention that 
■p p (]R 3 ) = {0}, if p < 0. Another relevant polynomial space is 

-Pp(divO,IR 3 ) := {qe T P (R 3 ) : divq = 0}. (3.1) 

Deep insights can be gained by regarding edge elements as discrete 1-forms. This pro- 
vides a very elegant construction of higher order edge element spaces and immediately 
reveals their relationships with standard Lagrangian finite elements and H(div, O)- 
conforming face elements (see below). In particular, the Poincare lifting becomes a 
powerful tool for building discrete differential forms of high polynomial degree. This is 
explored in [27,28], [29, Sect. 3.4], and [4, Sect. 1.4] in arbitrary dimension, using the 
calculus of differential forms. In this article we prefer to stick to the classical calculus 
of vector analysis, because we are only concerned with 3D. We hope, that, thus, the 
presentation will be more accessible to an audience of numerical analysts. Yet, the 
differential forms background has inspired our notations: integer superscripts label 
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spaces and operators related to differential forms. For instance, Wp(Ai) can be read 
as a space of discrete 1-forms. 

According to [27, Sect. 3], for any T G Ai, a G T, we can obtain the local space 

as 

WpCT) = 7> P (R 3 ) + R a (7> p (divO,M 3 )) . (3.2) 

Independence of a is discussed in [27, Sect. 3]. The representation (|3.2p can be estab- 
lished by dimensional arguments: from the formula (|2.2p for the Poincare lifting we 
immediately see that "P p (K 3 )-r-Ra('Pp(K 3 )) C W£(T). In addition, from [38, Lemma 4] 
and [27, Thm. 6, case ? = 1, n = 3] we learn that the dimensions of both spaces agree 
and are equal to 

dimV^CO = |(l+p)(3 + p)(4+p) . (3.3) 

As a consequence, the two finite dimensional spaces must agree. 

For the remainder of this section, which focuses on local spaces, we single out a 
tetrahedron T G M. . On T we can introduce a smoothed Poincare lifting according 
to (|2~3| with B = T and a suitable $ G Cfi°(T) complying with ([2^4)1 . An immediate 
consequence of (|3. 21) is that 

R T ('Pp(divO,M 3 )) c VVjJ(T) . (3.4) 

We introduce the notation !F m (T) for the set of all m-dimensional facets of T, m = 
0, 1, 2, 3. Hence, Tq(T) contains the vertices of T, J-\(T) the edges, ^(T) the faces, 
and Tz(T) = {T}. Moreover, for some F G T m {T), m = 1,2,3, V P (F) denotes the 
space of m-variate polynomials of total degree < p in a local coordinate system of the 
facet F, and T > P {F) will designate corresponding tangential polynomial vectorfields. 
Further, we write 

Wp(e) = Wp(T) ■ t e , *e the unit tangent vector of e, e € T\(T) , (3.5) 
W£(/) = Wp(T) x n f , n/ the unit normal vector of /, f£ T 2 {T) , (3.6) 

for the tangential traces of local edge element vectorfields onto edges and faces. Simple 
vector analytic manipulations permit us to deduce from (|3.2[) that 

W p 1 (e)=P p (e), ee^T), (3.7) 
K(f)=-P p (f) + R 2 a D (V p {f)), a€f, feT 2 (T), (3.8) 

where the projection R^ of the Poincare lifting in the plane reads 



R 2 a D (u)(x) := [ tu(a + t(x - a)](x - a) dt 
Jo 



ael 2 . (3.9) 



It satisfies div r R^ D (u) = u for all u G C°°(M. 2 ). We point out that, along with 
the formulas (|3.7[) and (|3.8p are special versions of the general representation formula 
for discrete 1-forms, see [27, Formula (16)]. Special facet tangential trace spaces will 
also be needed: 

Wp(e) :={ue W£(e) : Judl = 0}, e£^(T), (3.10) 

W p (/) := {u G Wlif) : u • n eJ = Ve G ^(T), e C 9/} , / G ^ 2 (T) , (3.11) 
Wj(T) := {u G Wp(T) : u x n/ = V/ G T 2 {T)} . (3.12) 
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Here rif represents an exterior face unit normal of T, n e ,f the in plane normal of a 
face w.r.t. an edge e C df. 

According to [38, Sect. 1.2] and [27, Sect. 4], the local degrees of freedom for 
Wp (T) are given by the first p — 2 vectorial moments on the cells of M. , the first p — 1 
vectorial moments of the tangential components on the faces of M. and the first p 
tangential moments along the edges of T, see (|3.14j) for concrete formulas. Then the 
set dofp(T) can be partitioned as 

dofi(T) = |J ldfi(e) U |J ldf p (/) U ldf p (T) , (3.13) 

where the functionals in ldf p (e), ldf p (/), and ldf p (T) are supported on an edge, face, 
and T, respectively, and read 

K G ldfp(e) => k(u) = f e p£- t e dl for e G Fi(T), suitable p G V p (e) , 

K G ldfp(/) =► k(u) = p • (£ x n) dS for / G .F 2 (T), suitable p e V p -i(f) , 
k e ldfp(T) =>■ k(u) := f T p ■ £dx for suitable p G V P - 2 {T) . 

(3.14) 

These functionals are unisolvent on W p (T) and locally fix the tangential trace of 
u G W£(T). There is a splitting of W£(T) dual to (|3"T5)) : Defining 

y p (F) := {v G W X (T) : k(v) = V K e dofJ(T) \ ldf*(F)} (3.15) 
for F G J- m (T), m — 1,2, 3, we find the direct sum decomposition 

w P x cn = £ E ( 3 - 16 ) 

m=l Fer m (T) 

In addition, note that the tangential trace of u G X^{F) vanishes on all facets ^ F, 
whose dimension is smaller or equal the dimension of F. By the unisolvence of dof p (T) , 
there are bijective linear extension operators 

EliP-W&^yfe), e G T\(T) , (3.17) 
E}, p :W;(/)-^(/), /e^ a (T). (3.18) 

The curl connects the edge element spaces W p (M) and the so-called face element 
spaces of discrete 2-forms [38, Sect. 1.3] 

W 2 p {M) := {v G H (div, Q) : v, T G W*(T) VT G X} , 
W p 2 (T) := {v G C°°(T) : v(x) = p(x) + ?(*)*, p G 7> P (T), g G 7> P (T)} . 

An alternative representation of the local face element space is [27, Formula (16) for 
I = 2, n = 3] 

W P 2 (T) = 7> P (T) + D„(P P (T)) , (3.19) 
where the appropriate version of the Poincare lifting reads 

(D a u)(x):=[ t 2 u(a + t(x - a))(x - a)dt , a E T . (3.20) 
Jo 
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Like p.2p this is a special incarnation of the general formula (16) in [27]. Again, dimen- 
sional arguments based on [38, Sect. 1.3] and [27, Thm. 6] confirm the representation 
[[3~20]) , We remark that div D a u = u, see [25, Prop. 1.2]. 
The normal trace space of W 2 (T) onto a face is 



W P 2 (/) := W P 2 (T) • n f = V p (f) , / G T 2 (T) 
and as relevant space "with zero trace" we are going to need 



(3.21) 



W 2 p (f) := {u G W p 2 (/) : J udS = 0} , / G F 2 (T) 
VV*(T) := {u G W 2 (T) : u • n 9T = 0} . 



(3.22) 
(3.23) 



The connection between the local spaces Wp(T), W 2 (T) and full polynomial 
spaces is established through a local discrete DeRham exact sequence: To elucidate 
the relationship between differential operators and various traces onto faces and edges, 
we also include those in the statement of the following theorem. There rif stands for 
an exterior face unit normal of T, n e j for the in plane normal of a face w.r.t. an edge 
e C df, and 4, is the differentiation w.r.t. arclength on an edge. 

Theorem 3.1. For f € ^{T), e e Fi(T), e C df, all the sequences in 



const 



const 



const 



P+i 



(T) 



- Pp+i(f) 

■4 



grad 



VVp(T) 



curl 



W P 2 (T) 



div 



w^/) VpU) 



{0} 



{0} 



{0} 



are exact and the diagram commutes. 

Proof. The assertion about the top exact sequence is an immediate consequence 
of representations (|3.2p and (|3.19[) and the relationships 



curlR a (u)=u Vu G -p p (div0,T), divD a (w)=u Vu G V P (T) . 

For further discussions and the proof of the other exact sequence properties see [27, 
Sect. 5 for n = 3]. 

4. Projection based interpolation. The degrees of freedom introduced above 
define local finite element projectors onto Wp(T). In conjunction with suitably defined 
interpolation operators for degree p Lagrangian finite elements, they possess a very 
desirable commuting diagram property [27, Thm. 13], which will be explained below. 
However, they do not enjoy favorable continuity properties with increasing p. Thus, 
L. Demkowicz [19-21], taking the cue from the theory of p- version Lagrangian finite 
elements, invented an alternative in the form of local projection based interpolation. 

4.1. Projections, liftings, and extensions. Again, consider a single tetra- 
hedron T G M. and fix the polynomial degree p € N. Following the developments 
of [29, Sect. 3.5], projection based interpolation requires building blocks in the form 
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of local orthogonal projections Pj, and liftings L^|. Some operators will depend on a 
regularity parameter < e < |, which is considered fixed below and will be specified 
in Sect. [5] To begin with, we define for every e S T\(T) 

P\ p : H-^(e) ^ |^+i(e) = VV*(e) (4.1) 

o 

as the ff -1+e (e)-orthogonal projection. Here, V P (F) denotes the space of degree p 
polynomials on a facet F that vanish on dF. 

Similarly, for every face / G J~2(T) introduce 

P} iP :H-^(f) » curl r ^ P+ i(/) = {v 6 VV p (/) : div r v = 0} , (4.2) 

Pl P :H-^(f) - divrW^(/) = W%f) , (4.3) 
as the corresponding if _3+e (/)-orthogonal projections. Eventually, let 

P^ p :L\T) i * gradP p+ i(T) = {v e Wj(T) : curlv = 0} , (4.4) 

P 2 Tp :L 2 (T) h-> curlW„(T) = {v e W 2 JT) : divv = 0} , (4.5) 



P 3 Tp :L 2 {T)^divW 2 JT) = {veV p (T) : [v(x)dx = 0} : (4.6) 

JT 



stand for the respective L 2 (T)-orthogonal projections. 
The lifting operators 



li p :W p (e)^V p+1 (e) , ee^i(T), (4.7) 

n n 

4, 



{v G Wj(/) : divrv = 0} » V p+1 (f) , f e F 2 (T) , (4.8) 
{v 6 W p (T) : curlv = 0} -> £ p+1 (T) , (4.9) 



are uniquely defined by requiring 

j^lpu^u VueWj(e), (4.10) 

curl r L} p u = u Vue {W p (/) : div r v = 0} , (4.11) 
grad UppU = u Vu e {v e W p (T) : curlv = 0} . (4.12) 

Another class of liftings provides right inverses for curl and divr: Pick a face / 6 
•^(T), and, without loss of generality, assume the vertex opposite to the edge e to 
coincide with 0. Then define 

,2 . f div r W p (/) ^ WUf) , , 

f ' p -\ « - Ro u-curlrES-^^u.^). (4 ' 13) 

This is a valid definition, since, by virtue of definition (|3.9p , the normal components 
of Rq D u will vanish on df \ e. Moreover, divrRo D u = u ensures that the normal 
component of Rq D u has zero average on e. We infer 

(curl r E| ip Li p ((R^u • t» ?i/ ),-) • n ?l/ ) = 

|l-l, P ((Ro D u) • »e,/) - Rg D u • ne,/ on e , 



4 The parameter £ in the notations for the extension operators E^, the projections P"., and the 
liftings l_'„ refers to the degree of the discrete differential form they operate on. This is explained 
more clearly in [29, Sect. 3.5]. 
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and see that the zero trace condition on df is satisfied. The same idea underlies the 
definition of 

2 f curlW p (T) i ^ W p (T) 
L t,p : | u ~ R u-gradE^LL p (((R 0U )xn 7 ) |7 ), ( 4 ' 14 ) 

where / is the face opposite to vertex 0, and the definition of 

J div W»(T) _> WJ(T) 
L t,p-| u ^ Don- curl EL j) L| ip ((D 0U -n 7 ) |J ). ( 4 - 15 ) 

The relationships between the various facet function spaces with vanishing traces can 
be summarized in the following exact sequences: 

{0} -J*-> V P+1 (T) -f=^ W p (T) -A Wl(T) -fU p p (T) -J*-, {0}, 

T,p T,p T,p 

{0} -^U wl(/) ^ {0}, 



L 2 



{0} V p+1 (e) -f-^ P p (e) {0}, 



(4.16) 



where V P {F) designates degree p polynomial spaces on F with vanishing mean. These 
relationships and the lifting mappings are studied in [29, Sect. 3.4]. 
Finally we need polynomial extension operators 

E° e , p :V p+1 (e)~V p+1 (T), (4.17) 

E° Lp :r p+1 (f)^r p+1 (T) (4.18) 

that satisfy 

E° pV =0 Ve'e^i(T)\{e}, (4.19) 

E?, pV =0 V/' e ^ 2 (T) \ {/} . (4.20) 

Such extension operators can be constructed relying on a representation of a polyno- 
mial on F, F <G T m {T), m — 1,2, as a homogeneous polynomial in the barycentric 
coordinates of F, see [29, Lemma 3.4]. As an alternative, one may use the polynomial 
preserving extension operators proposed in [22,37] and [1]. We stress that continuity 
properties of the extensions E' F , I = 0, 1, F S T m {T), are immaterial. 

4.2. Interpolation operators. Now, we are in a position to define the projec- 
tion based interpolation operators locally on a generic tetrahedron T with vertices Oj, 
i = 1,2,3,4. 

First, we devise a suitable projection (depending on the regularity parameter 
< e < |, which is usually suppressed to keep notations manageable) 

nU= n T») : C°°(T) ~ P p+1 (T) (4.21) 
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for degree p Lagrangian i/ 1 (il)-conforming finite elements. For u E C°(T) define (Aj 
is the barycentric coordinate function belonging to vertex of T) 



«(°) 



u-^u{ ai )\, (4.22) 



i=l 



:=u>(°) 

««:=«(»)- £ E° p L^P^^(°) |e , (4.23) 

S v ' 

:=iu(i) 

U ( 2 ) : = U «- ]T E UL P /, P <™M« (1) |/) . (4-24) 

/e^i(T) 



II^pU := L^ p P^ pg radu (2) + w (2) + + w (0) . (4.25) 

Observe that w^\ F = for all F e F m (T), < m < i < 3. We point out that u/ ) is 
the standard linear interpolant of u. 

Lemma 4.1. 77ie linear mapping n^ p , p e No, is a projection onto Cp p+ \(T) 

Proof. Assume u £ V p +i(T), which will carry over to all intermediate functions. 
Since u^°\zi) = 0, i = 1, . . . , 4, we conclude from the projection property of P* that 
LgP;!,^;u( )| e = u^°'| e for any edge e € ^(T 1 ). As a consequence 



^ E^u^ie =* «[?=(> Vee^i(T). (4.26) 

ee^i(T) 



We infer Lj. p Py.curlr(M < - 1 ^/) = «^|/ on each face / G ^(T), which implies 



EU uW \f) =* « (2) l/ = V/e^ 2 (T). (4.27) 



This means that p P^ p gradt^ 2 ' = and finishes the proof. □ 

A similar stage by stage construction applies to edge elements and gives a pro- 
jection 



n^(= n^(e)) : C°°(T) W^T) : (4.28) 
for a directed edge e := [a^a,-] we introduce the Whitney- 1-form basis function 

b e = A 4 grad Xj - Xj grad A l . (4.29) 
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These functions span W^T). Next, for u E C°(T) define 

u<°> := u - ( / u • ds) b e , (4.30) 



U (D := u (0) 



u (2) . = U (D 



u (3) . = u (2) 



:=w<°) 

]T gradE° e)p L^P^((u(°)-t e ) |e ), (4.31) 

~ w (i) 

£ E} ip L 2 /ip P^div r ((u( 1 ) x , (4.32) 

:=w< 2 > 

2 gradE°. p L} !?) P} !p ((u( 2 ) x n,),,) , (4.33) 

/6^2(T) 



:=w<3) 

u (4) . = u (3) _ |2 p2 cur j u (3) ^ ( 4 34 ) 



II^u := grad L^P^u^ + w' 4 ' + w< 3 > + w< 2 > + w« + w' ' . (4.35) 

The contribution w(°) is the standard interpolant 11^ of u onto the local space of 
Whitney-l-forms (lowest order edge elements, see [35, Sect. 5.5.1]). The extension 
operators were chosen in a way that guarantees that w^ 2 ) ■ t e = and w^ 3 ) • t e = 
for all e G Fi(T). 

Lemma 4.2. The linear mapping Ib^ p , p € No, is a projection onto W p (T) and 
satisfies the commuting diagram property 

II^ p o grad = grad oU^ p on C°°(T) . (4.36) 



Proof. The proof of the projection property runs parallel to that of Lemma 14.11 
Assuming u S W p (T) it is obvious that the same will hold for all uW and wW from 
(|4.30[) - (|4.35[) . In order to confirm that all projections can be discarded, we have to 
check that their arguments satisfy conditions of zero trace on the facet boundaries 
and, in some cases, belong to the kernel of differential operators. 

First, recalling the properties of the interpolation operator IIq for Whitney-l- 
forms, we find (u<°) • t e ) |e e W*(e). This implies 

grad E° p L^ P U ( U(0) ' *e)|e) = ( U(0) ' *e)|e Ve G ^( T ) . ( 4 ' 37 ) 



and 



(u«.t o ), e = Vee^i(T) 



(4.38) 
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We see that (uW x n f ), f G for any / G T 2 {T), so that 

P^div r ((u« x n f ) lf ) = div r ((u (1 ) x (4.39) 

divrL2 p P 2 /p divr((u( 1 ) x n/),,) = dW r ((u^ x n/)^) (4.40) 

divr((u( 2 ) x nf ) ]f ) = V/ G T 2 {T) , (u< 2 > • t e ) |e = Ve G ^i(T) (4.41) 

P}, p ((u( 2 ) x n f ) lf ) = (uW x n,),, V/ G :F a (T) (4.42) 
=> grad E° /!p L} !p P} !p ((u( 2 ) x n,),,) x n f = (u< 2 ) x n,),, V/ G T 2 (T) (4.43) 

=► (u< 3 > x n f ) lf = V/ G F 2 (T) (4.44) 

=> P| p curl u (3) = curl u (3) (4.45) 

=> curl Lf, p Pf, p curl u (3) = curl u (3) (4.46) 

=> curlu (4 >=0 =► Pj, U W=uW (4.47) 

=► gradL^P^uW =u< 4 ) , (4.48) 

which confirms the projector property. 



Now assume u = gradu for some u G C°°{T). The commuting diagram property 
will follow, if we manage to show gradu' ) = u^°\ gradu' 1 ' = u' 1 ), gradu' 2 ) = u' 3 \ 
etc., for the intermediate functions in (|4.22| - (|4.25| and (|4.30|) - (|4.35|) . respectively. 

By the commuting diagram property for the standard local interpolation operators 
onto the spaces of Whitney-O-forms (linear polynomials) and Whitney- 1-forms, we 



conclude 

gradu' ' = u(°) => ^« (0) |e = (u (0) • t e )\ e Ve G ^i(T) (4.49) 

=> u (1) = gradu (1) => div r ((u (1) x n f ), f ) = V/GJ" 2 (T) (4.50) 

u( J )=u« (4.51) 

=> (u (2) x n f ) lf = curl r u (1) / V/ G T 2 {T) => u (3) = grad u (2) (4.52) 

uW = u < 3 ) . (4.53) 



Of course, analogous relationships for the functions ioW and hold, which yields 
n^ p u = gradllPr pU. □ 

Following [29, Sect. 3.5], a projection based interpolation onto W 2 (T), the oper- 
ator n| p (= n| >p (e)) : C°°(T) i-> W 2 (T), involves the stages 

u(°' :=u-( [ u - n f dS ) h f, ( 4 - 54 ) 



u W:=uC°J- £ curlE /,P L L P L(( u(0) ■«/),/) (4-55) 



u (2) . = u (l)_ L 3 p p3 pdivu (l) (4 56) 

" : ^ ' 

:=w< 2 > 

n^ p u := curl L^P^u^ + W ( 0) + w« + w< 2) . (4.57) 
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Here, refers to the local basis functions for Whitney-2-forms [29, Sect. 3.2]: 

bf = Xi grad Xj x grad Xk + Xj grad Xk x A, + Xk grad Aj x Xj . (4.58) 
Analogous to Lemma T4. 2 1 one proves the following result. 

Lemma 4.3. The linear operator 11^ p , p e No, is a projection onto W p (T) and 
satisfies the commuting diagram property 

Ii 2 T p o curl = curl oII^ p on C°° (T) . (4.59) 

The next lemma makes it possible to patch together the local projection based 
interpolation operator to obtain global interpolation operators 

nj, : C°°{T1) ^ W l p (M) , 1 = 1,2. (4.60) 

Lemma 4.4. For any F G T m (T), m = 0,1,2, and u G C°°(T) the restriction 
H^ p u^ F depends only on u\p. 

For any F G J- m (T), m — 1,2, and u G C°°(T) the tangential trace of Yl\ p \i 
onto F depends only on the tangential trace of u on F. 

For any face f G ^(T) andu G C°°(T) the normal trace ofUj, u onto f depends 
only on the normal component of u on f . 

Proof. The assertion is immediate from the construction, in particular, the prop- 
erties of the extension operators used therein. □ 

It goes without saying that density arguments permit us to extend lip, I = 0, 1, 2, 
to Sobolev spaces, as long as they are continuous in the respective norms. (Repeated) 
application of trace theorems [26, Sect. 1.5] reveals that it is possible to obtain con- 
tinuous projectors 

U° p : H 1+S (fl) ^{ve H\n) : v\ T £ V P+1 (T) VTeM}, (4.61) 

Il£ : H^ +S (fl) ^ W$(M) , (4.62) 

U 2 p : H S (Q) ^>W 2 p {M) , (4.63) 

for any s > 1. In addition, by virtue of Lemma l4.41 zero pointwise/tangential/normal 
trace on dQ of the argument function will be preserved by 11^, I = 0, 1, 2, for instance, 

nJ(H' +s (fi)nH (curl,n)) = W*(M)nH (curl,n) . (4.64) 

5. Interpolation error estimates. Closely following the ingenious approach 
in [21, Section 6] we first examine the interpolation error for LTy. . Please notice that 
IIj, still depends on the fixed regularity parameter < e < | . The argument function 
of II^p is assumed to lie in H 1+S (T) for some s > |, cf. (|4.61|) . The continuous 
embedding H 1+S (T) <^-> C°(T) plus trace theorems for Sobolev spaces render all 
operators well defined in this case. 

We start with an observation related to the local best approximation properties 
of the projection based interpolant. 

LEMMA 5.1. For any u G H 1+S (T) holds 

(grad(u-n° :P u),grad«) i2(T) =0 Vw G V P +i(T) , (5.1) 
(cuvl r (u-U° Tp u) curl r v) =0 W G V P+ i(f), f G F 2 (T) , (5.2) 

\ U / H 2 (/) 

(|( u - n r, P «)| e ,|«J ffi+ =0 VveV p+1 (e), eG^(T). (5.3) 
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Proof. We use the notations of (|3~2"2"|) - (pl~2l))) . Setting w := w (0) + + w (2) , we 

find 

II^u = L^ iP P^ p grad(u - w) + w , (5.4) 
which implies, because is a right inverse of ^ , 

grad Ily pii = grad it + (Id — Pr )P ) grad w . (5.5) 

This means that m— gradn^ p u belongs to the range of Id— P^ p and (|5.ip follows from 
(|4.4p and the properties of orthogonal projections. Similar manipulations establish 
JO): 



curlrn^ jp u^. = curlr?«|/ (5.6) 
= curl r L} p P}^curl r u (1) + curl r (w (0) + w^)\ f (5.7) 



= ld 

= P^curlrtti/ + (Id - P} iP )curl r (^ (0) + u> (1) ) V/ G ^ 2 (T) . (5.8) 

The same arguments as above verify (I5.3[) . □ 

From this we can conclude the result of [21, Section 6, Corollary 1]. To state it 
we now assume a dependence 

° <C = e(p):= 101ogUl) < ^ ^ 

of the parameter e in the definition of the local projection based interpolation opera- 
tors. Below, all parameters e are linked to p via (|5.9[) . Please note that we retain the 
notation (n^) gls} , I = 0, 1, 2, for these new families of operators. 

Theorem 5.2 (Spectral interpolation error estimate for 11^ ). There is a con- 
stant Ct > depending only on T and \ < s < 1, and, in particular, independent of 
p, such that 

\(\d-n^ p )v\ Hl{T) <C T l -^P\v\ Hl+s{T) Vu e H 1+S (T) , „>1. (5.10) 

Stable polynomial extensions are instrumental for the proof, which will be post- 
poned until PageHH First, we recall the results of [37, Thm. 1] and [1, Thm. 1]: 

Theorem 5.3 (Stable polynomial extension for tetrahedra). For a tetrahedron T 
there is linear operator St ■ Hz{dT) i— > H l (T) such that 

S T u ldT = u VueH*(dT), (5.11) 
\S T u\ H1{T) <C\u\ Hi{dT) VueHHdT), (5.12) 
S T weV p+1 {T) VweV p+1 (T) ldT , (5.13) 

where C > only depends on the shape regularity measur^ of T . 

Theorem 5.4 (Stable polynomial extension for triangles). Given a triangle F, 



5 The shape regularity measure of a tetrahedron is the ratio of the radii of its circumscribed sphere 
and the largest inscribed sphere. 
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there is a continuous linear mapping S F : L 2 (dF) i-> H?(T) such that 

\S F u\ H1(F) <C\u\ Hi(gF) Vu€Hi(dF), (5.14) 
S F weV p+1 (F) VweV p+i (F) lgF , (5.15) 

where C > depends only on the shape regularity measure ofT. 

By interpolation in Sobolev scale from the last theorem we can conclude 

3O0: \S F u\ HS{F) <C\u\ HS _i (dF) VueH s -HdF), i < s < 1 . (5.16) 

We also need to deal with the awkward property of the _ff ^ (<9T)-norm that it 
cannot be localized to faces. To that end we resort to a result from [34, Proof of 
Lemma 3.31], see also [21, Lemma 13]. 

Lemma 5.5 (Splitting of H ~s (<9T)-norm). The exists C > depending only on 
the shape regularity of the tetrahedron T such that 

\u\ i +e <- V [til Vu e H^ +e (dT), < e < - . (5.17) 

/£;f 2 (t) 

Another natural ingredient for the proof are polynomial best approximation esti- 
mates, see [40] or [37, Sect. 3]. 

Lemma 5.6. Let F be either a tetrahedron or a triangle. Then, there is a constant 
C > depending only on F such that for all p > 1 

inf \u-v p \ Hr(F) <Cp r - 1 - s \u\ H1+ . {F) Vu G H 1+8 (F), < r < 1 . (5.18) 

Define a semi-norm projection Qt, p ■ -ff 1 (T) i— > V p -\-i(T) on the tetrahedron T by 

/ grad(w - Q T , P u) ■ gradu p da; = Vw p G V p+ i(T) , 
Jl (5.19) 

it — Qt, p u dx = , 

and semi-norm projections Q/, p : H s+ ^(f) i— > V P +i(f), f G ^{T), by 
(curl r (tt- Q/, P tt),curl r i;p) ffe _i =0 Vw p G , 



y it — Qj ^udcc = 



(5.20) 



These definitions involve best approximation properties of Qt, p u and Qf iP u. Thus, we 
learn from Lemma I5.6I that with constants independent of0<e<i<s<l 

\u-Q TiP u\ H1(T) <C(p+l)- s \u\ H1+s{T} VueH s (T), (5.21) 
\u - Q f , p u\ Hi+t(f) < C(p + iy- s \u\ Hi+S{T) Vu G H?+ S (f) . (5.22) 

The latter estimate follows from the fact that |-| i .,, and ||curlr-|| ii,, . are 
equivalent semi-norms, uniformly in e. 
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We also need error estimates for the L 2 (e)-orthogonal projections, 

Ql p :L 2 (e)^V p+1 (e) , e E T\{T) . (5.23) 



Lemma 5.7 (see [21, Lemma 18]). With a constant C > independent of p, 
< e < \, and 2e < r < 1 + e 

e-Ql p u\ Hc[e) <C(p + l) 2£ ->|^ (e) VuGffCeJnffoHe) . 



Proof. Write l ejP : i?o(e) i— > Pp+i for the interpolation operator 

(U, P «)(0=«(0) + y o (Q e ,p^)(r)dr, 0<C<|e|, 

where £ is the arclength parameter for the edge e and Q eiP : L 2 (S1) i— > V p (e) is the 
L 2 (e)-orthogonal projection. From [40, Sect. 3.3.1, Thm. 3.17] we learn that 

\u-\ e ,pu\ a i(e)^ C <P+ 1 )~ 1 \ u \H'W VueH 2 (e), (5.24) 
\\u-\e, P u\\ L3(e) <C(p+l)- m \u\ Hm{e) VueH m (e), m=l,2. (5.25) 

Here and the in the rejnainder of the proof, all constants may depend only on the 
length of e. As \ e<p u £ V p +i(e) for u € Hg(e), this permits us to conclude 

\\u-Ql p u\\ L2(e) < \\u-\ e ,pu\\ LHe) <C(p+l)- l \\u\\ BHe) , (5.26) 

which yields, by interpolation between H (e) and L 2 (e), 

\\u-Ql p u\\ L2(e) <C(p + l)-"\\u\\ Hq(e) , 0<<?<1, (5.27) 

where C > is independent of q. On the other hand, using the inverse inequality [6, 
Lemma 1] 

Hl ffl(e) <C(p+l) 2 || U || i2(e) V«eP p+1 (e) (5.28) 
and (I5.24p . (|5.25|) we find the estimate 



< |« - l e , p u| ff i (e) + (P + I) 2 ||Qe,p« - kpHL^e) 

< \U - \ e ,pU\ H1{e) + C(p + l) 2 \\u - l^pWH i2(e) 
<C|HI^(e). 

r-2e 



(5.29) 



Interpolation between (I5.27|) with q = ^rzr an d (j5 . 29[) finishes the proof. □ 

Proof, [of Thm. 15.21 borrowed from [21, Sect. 6]] Orthogonality (15. ip of Lemma l5.1l 
combined with the definition of Qt, p involves 

J grad((n^ - Qt, p )u) • grad Up da; = Vv p £ V P +i(T) . (5.30) 

T 
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Hence, (n!^ p — Qx, p )u turns out to be the |-|^i ^-minimal degree p + 1 polynomial 
extension of (IPr p — Qt,p)u^ t , which, thanks to Thm. [5~31 implies 



(n° TjP - Q T , p )u\ H < S T ((U° T!p u - Qt, p u) ) 



ff!(T) 



< c 



(U° Tp u - Q T , P u) 



\dT 



H?(dT) 



Thus, by the continuity of the trace operator i? 1 (T) i— > H2 (dT), 
\u-U° Tp u\ HHT) < \u-Q T , p u\ HHT) 



C 



DT 



H2(dT) 

< C ( \u - Q T , P u\ H1(T) + (u - U^ p u) 



(u - Qt, p u),q T 



(5.31) 



H2(dT) 



P >\dT 



(dT) 



(5.32) 



To estimate 



(«- U T,p u ) ldT 

(u - n° u ) 



H2~ (dT) 



we appeal to Lemma 15751 and get 



\0T 



< 



H2(dT) 



(u - n° T!P u) 



\8T 



H2 +e (dT) 



C 



<- J2 k w - n V) 



1/ 



(5.33) 



H1+'(f) 



Next, we use (|5.2j) from Lemma 15.11 together with (|5.20j) . which confirms that 
(Tl T v u), — Qf, P u is the minimum |-| i + -seminorm polynomial extension of 

(Ll^ p u)^ d j — Qf^ p (u)\ g f. Hence, based on arguments parallel to the derivation of (|5.32[> . 
this time using Thm. 15.41 we can bound 



(u - U T , p u) lf K i+ < hi/ - Q/,P u | H i+. (/) + C ( u t, p u ~ Qf, P u) 



I Of 



H*(df) 

(5.34) 



where the (e-independent !) continuity constant of the trace mapping S/ enters the 
constant C > 0. Also recall the continuity of the trace mapping H^ +e (f) i— > H e (df) 
[34, Proof of Lemma 3.35]: with C > independent of e, 



C 



h8/L.( 8/ )<-^NI ff i + . (/) Vueff> +e (/) 



(5.35) 



Use this to continue the estimate ()5.34[) 

1 



(u-U^ p u) ]f 



H? +e (f) 



< c 



-j=\u lf -Q Lp u\ H i 



2+£ (/) 



( u ~ U T,pU) ldf 



H'(df) 



(5.36) 



As e < ^, we can localize the norm 



(«-n^) |a/ 



H'(df) 



to the edges of /: 



( u - U T,pU) ldf 



< 



C 



eeJFi(T),eC9/ 



H«(e) 



(5.37) 
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Recall the e-uniform equivalence of the norms |-|jje( e ) and || ;|§"||#-i+ e ( e y Hence, owing 
to (|5.3j) . we have from Lemma 15 . 71 with r = s: 



H«(e) 



< C inf 



(u~U%u) ~v p 



H'(e 



< C 



(u-n%u) le -Ql p ((u-IL%u) ]e ) 



H'(e) 



(5.38) 



<c( P + i) 2e - s (u-n° T0 u). 



Moreover, H 1+S (T) is continuously embedded into C°(T). Consequently, applying 
trace theorems twice and appealing to the equivalence of all norms on the finite 
dimensional space V\ (T) , 



(u — Tlx u), 



< \ u \e\ it,, , 



HHe 



— ^\ u \h 1 + s (T) i 



(5.39) 



where C > depends on s and T, but not on p. Combining the estimates (|5.32|) . 
(15331) . (j5~36)l . and (j5~37)l . (Q8ll with (j5~39)l . we find 

| U - n T,p U Ui(T) " C (\ U -QT,pU\ H1(T) + -^ \ U \f-QfA u \f)\ H i+c u) + 

( P + l) 2t - s ^ ^ 
e(i-e) ^ 



76^2 (T) 



lffi+=(T) J i 



(5.40) 

with C > independent of p. Finally, we plug in the projection error estimates (|5.21D . 
(|5.22|) . and arrive at (C > independent of u, e, p, s) 



\u - U° T Je)u\ Hl(T) < C( (p + l)- s \u\ H1+s(T) + 



(p+1)- 



,3/2 



/e^2(T) 



'(/) 



v2 ee^i(T) 

(5.41) 

The choice (|5.9|) of e together with an application of trace theorems then finishes the 
proof. □ 

The next lemma plays the role of [8, Lemma 9] and makes it possible to adapt 
the approach of [8, Sect. 4.4] to 3D edge elements. 

Lemma 5.8. // | < s < 1 and u S H S (Q) satisfies curlu T S V P {T) for all 
T e M, then 



|(id-n>| 



L 2 (0) 



< c 



log 3 / 2 



(H\ w 



(O) 



curl u| 



L 2 (T)i 



(5.42) 



iwt/i a constant C > depending only on M. and s. 

Proof. Pick any u complying with the assumptions of the lemma. The locality of 
the projector allows purely local considerations. Single out one tetrahedron T E A4, 
still write u = U| T , and split on T 

u = (u — Rt curl u) + Rt curl u , (5.43) 

Note that the properties of the smoothed Poincare lifting Rt stated in Thm. l2~3l imply 
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1. curl(u — Rx curlu) = on T, as a consequence of (I2.8p . and 

2. R T curlu G ff^T) and the bound 

||R T curlu|| ffl(T) < C||curlu|| i2(n) , (5.44) 

where here and below no constant may depend on u or p. 
Hence, as u G H S (T), there exists w G H 1+S (T) such that 

u = grad u + Rt curl u . (5.45) 

The continuity of Rt reveals that, with a constant C > depending only on T, 

\v\ H i+s {T) < ||u|| H3(T) + |R T curlu| ffl(T) < ||u|| Hs(T) +C*||curlu|| L2(T) . (5.46) 

By the assumptions of the lemma and (|3.4[) we know that 

R r curl u G W$(T) . (5.47) 

By the commuting diagram property from Lemma 14.21 and the projector property of 
n^ p the task is reduced to an interpolation estimate for n!^ p : 

(Id - n^)u grad(ld - n° Tp )v + (Id - n^ p )R T curlu (5.48) 

v -v ' 

=o 

As a consequence, invoking Thm. 15.21 



||( ld - n T :P ) u || L 2 (r) \O d - U T,p)v\ H1{T) 

<C-^\v\ m+s{T) < C-^(||u|| HS(r) + ||curlu|| i2(T) ), (5.49) 

which furnishes a local version of the estimate. Squaring (|5.49|) and summing over all 
tetrahedra finishes the proof. □ 

6. Discrete compactness. Smoothness of the solenoidal part of the Hclmholtz 
decomposition of if (curl, f2) and ifo(curl, Q), respectively, plays a pivotal role. It 
can be deduced from elliptic lifting theorems for 2nd-order elliptic boundary value 
problems [18, Ch. 6]. Proofs of the following lemma can be found in [29, Sect. 4.1] 
and [2, Sect. 3]. 

Lemma 6.1. For any Lipschitz polyhedron C M 3 there is a | < s < 1 such that 
X := ffo(div, fi) n if (curl, f2) and X :— if(div, O) n ifo(curl, f2) are continuously 
embedded into ff s (f2), that is, 

3C = C( Si n)>0: ||u|| H . (n) <c(||u|| ff(curiitl) + ||u|| K(dlVin) ) Vu G X . (6.1) 



We first verify the discrete compactness property of Def. Il.ll for e = 1: consider a 
sequence (u p ) pe N, which satisfies 

(i) u p G Wl(M) , (6.2) 

(ii) (u p ,z p ) L2{n) =0 Vz p G {v G W$(M) : curlv = 0} , (6.3) 

(iii) KIIl^o) + ||curlu p || L2(n) < 1 Vp G N . (6.4) 
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Theorem 6.2. A sequence (u p ) pe ^ compliant with (|6.2p - (|6.4p possesses a subse- 
quence that converges in L 2 (Q). 

Proof. The proof resorts to the "standard policy" for tackling the problem of 
discrete compactness, introduced by Kikuchi [32,33] for analyzing the h- version of edge 
elements. It forms the core of most papers tackling the issue of discrete compactness, 
see [9, Thm. 2], [8, Thm. 11], [29, Thm. 4.9], [23, Thm. 2], [8, Thm. 11], etc. 

We start with the continuous Helmholtz decomposition of u p : let u p be the unique 
vector field in if (curl, O) with 

curl Up = curl u p , (6-5) 
(u p , z) L2(0) =0 Vz G Ker(curl) n H (curl, Q) . (6.6) 

The inclusion grad_ff 1 (fi) C Ker(curl) enforces 

divUp = in , u p ■ n = on dfl . (6-7) 
Hence, by virtue of Lemma 16. 11 u p satisfies 

u p EH s (n) , \\u p \\ H3{n] <C\\u p \\ Hicurl ^ , (6.8) 

where C > depends only on and ~ < s < 1. 

In addition, u p is L 2 (f2)-orthogonal to Ker(curl) n ii(curl, fi), see (|6.6I) . Thus, 
using Nedelec's trick [38] , we have 



u pIIl 2 (o) — ( u p u pi u p npup + npup u p) £ 2( f2 - ) 

= (Up — Up, Up — n p Up) L 2/m ! 



(6.9) 



because, thanks to Lemma 1431 and (|6.5|) . 

curl(n p u p - Up) = (n 2 - Id) curl Up = . (6.10) 

=S> np5p-Up€{ve WjJ(A4): curlv^O}. (6.11) 
Hence, appealing to Lemma f5T5l with C > independent of p, 



I u p - u p\\mn) < K ~ K u p\\mn) ^ C ps (1 u p1h»(o) + llcurl u P || i2(n) ) 



log 3/2 p 

< C - s Kllff(curl,n) -> fOT P -» 00 ' 



(6.12) 



Since bounded in Ji s (Sl), by Rellich's theorem (u p ) pgN has a convergent subsequence 
in L 2 ($}): owing to (|6.12[) . the same subsequence of (up) pgN will converge in L 2 (Q). □ 

Theorem 6.3. Replacing W^M) withW£(M)nH (curl,Q,) in . the 

assertion of Thm. \6.2\ remains true. 

Proof. Since the projection based interpolation operators respect homogeneous 
Dirichlct boundary conditions, cf. (|4.64|) . the proof runs parallel to that of Thm. ROl 
We point out that now u p x n — on dfl instead of u p • n = 0, but Lemma T6. II can 
still be applied. □ 
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Now we are able to switch from e = 1 to general dielectric tensor, thus completing 
the proof of the main theorem Thm. 11.21 

Proof, [of Thm. O in the Introduction] We adapt the proof of [29, Thm. 4.9]. 
Consider a Jf(curl, £!)-bounded sequence (w p ) pgN , w p e Wp(M), in the Ll(fl)- 
orthogonal complement X^{M.) (see (11.21) ) of the discrete kernel of curl, i.e., 

(ew p ,z p ) i2(sl) =0 Vz p e Ker(curl)n W$(M) . (6.13) 

We continue with the L 2 (f2)-orthogonal discrete Helmholtz decomposition 

Wp = Up® L 2w£, u„ 6 W^jW), w° 6 Ker(curl) n Wy 1 ^) ■ (6.14) 

As (u p ,z p ) L2 ^ — for all z p e Ker(curl) n Wp(M), by Thm. 16.21 we can find a 
subsequence, again denoted by (u p ) pgN , with 

u p -2=^ q i n £ 2(q) _ ( 6 .i5) 

Since w p satisfies (I5.4ip . we conclude 

(ew°,Zp) L2(n) --(eu p ,zp) £2(n) Vz p 6 Ker(curl) n Wp(M) . (6.16) 

This can be regarded as a perturbed spectral Galerkin approximation of the following 
continuous variational problem: seek y € i?(curl0,il) := Ker(curl) n ff(curl, fi) 
such that 

( £ y> z ).L=(n) = - ( e{ h z )t'(n) Vz e fr(curlO,ft) , (6.17) 

which, obviously, has unique solution y. From Strang's lemma [14, Thm. 4.4.1] we 
infer 

||y-w°L (n) < c ( Zpew , ( ^Kc r( cu r ,) l|y - Zpll ^w + K-qii.W . ( 6 - 18 ) 

— >0 for p — >oo 

where C > depends only on e. Next recall, that there is a representation 

y = gradw + h, v e H\Q) , h G ?i 1 (f2) , (6.19) 

with H (fi) standing for the finite dimensional first co-homology space W 1 ^) of har- 
monic vector fields, which is contained in i?(curl 0, O)ni?o(div 0, ft), [29, Lemma 2.2] 
and [2, Prop. 3.14, Prop. 3.18]. Owing to Lemma [6.11 it belongs to H s (£l) for some 
s > i, which implies, thanks to Lemma [ 



h - riph|| i2(0) < C^^|h| H3(fi) - I) fc- - x 10.20) 



Further, the commuting diagram property of Lemma 14.31 confirms that curlllph = 0. 
Besides, asymptotic density of the spectral family of Lagrangian finite element spaces 
in i? 1 (fi) means that also the first term on the right hand side of (|6.18[) tends to zero 
as p — > oo . 

Thus, selecting the same subsequence of (w p ) p( _ N (and keeping the notation), it 
is immediate that 

Wp q + y ini 2 (f>). (6.21) 

The case of w p £ Wp(M) fl Ho(curl,il) is amenable to almost the same proof: the 
boundary conditions are imposed on all fields and in the counterpart of (|6.19|) the 
second co-homology space 7i 2 ($V) has to be considered [29, Lemma 2.2]. □ 
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